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ABSTRACTS  OF  PRESENTATIONS 


Molecular  Dynamics  Simulations  of  the  Phase  Transitions  of  the  Homopolymer 
of  p-hydroxybenzoic  Acid 

Michael  Nolan 

Materials  Modelling  Group,  National  Microelectronics  Research  Centre,  UCC,  Cork. 

We  present  the  results  of  force-field  calculations  and  molecular  dynamics  simulations  on  the 
homopolymer  of  the  liquid  crystal  polymer  p-hydroxybenzoic  acid.  A  model  unit  cell  with 
two  chains  and  six  monomers  per  chain  was  constructed  based  on  experimental  results. 
The  FF  calculations  were  used  to  optimise  the  lattice  structure.  The  MD  simulations  were 
carried  out  within  the  NPT  ensemble,  allowing  for  changes  in  the  cell  parameters  and  shape. 
Discontinuities  in  the  isobaric  specific  heat  capacity  as  a  function  of  temperature  were 
used  as  the  signature  for  the  phase  transitions.  Good  agreement  with  the  experimentally 
quoted  phase  transitions  was  found.  The  structural  changes  which  occur  at  the  first  phase 
transition  were  also  modelled.  We  also  present  an  analysis  of  the  molecular  motion  which 
occurs  at  and  determines  each  phase  transition. 


Molecular  simulation  of  the  separation  of  gas  mixtures  using  thin  nanoporous 
membranes 

Don  MacElroy 

Department  of  Chemical  Engineering,  NUI,  Dublin. 

A  nonequilibrium  molecular  dynamics  (NEMD)  simulation  technique  is  employed  to  inves¬ 
tigate  the  transport  of  binary  mixtures  of  hydrogen  and  methane  through  a  model  carbon 
membrane  of  varying  thickness.  Both  forced  flow  simulations  under  a  pressure  gradient 
and  isobaric  counterdiffusion  simulations  are  described  in  this  presentation.  The  principal 
conclusions  of  these  studies  are:  (i)  pore  entrance/exit  effects  may  need  to  be  taken  into 
consideration  in  membrane  design  for  hydrogen/hydrocarbon  separations  if  the  length  of 
the  controlling  pores  within  the  carbon  membrane  is  of  the  order  of  one  tenth  of  a  micron 
or  less;  (ii)  viscous  (convective)  flow  contributions  to  the  fluxes  of  the  individual  compo¬ 
nents  of  the  mixture  should  be  negligible  for  the  carbon  membranes  currently  in  use;  and 
(iii)  the  cross-coefficients  of  diffusion  appear  to  play  a  relatively  minor  role  in  the  normal 
(pressure  driven)  hydrogen/methane  membrane  separation  process  however  under  isobaric 
conditions  the  simulation  results  suggest  that  hydrogen/methane  cross-coupling  plays  a 
significant  role  in  hindering  the  hydrogen  counterdiffusion  flux  within  carbon  membranes 
containing  long,  narrow  pores. 


3 


An  investigation  into  surface  ordering  in  liquid  n-alkanes  with  molecular  dy¬ 
namics 

Paul  Smith 

Atomistic  Simulation  Group,  Queen’s  University,  Belfast. 

The  unusual  surface-ordered  state  of  liquid  n-alkanes  was  investigated  by  simulations  of 
a  slab  of  liquid  heptadecane  with  two  free  surfaces.  The  results  confirm  that  the  surface- 
ordered  state  consists  of  a  monolayer  of  molecules  aligned  perpendicular  to  the  surface 
which  are  mostly  in  an  all  trans  conformation.  The  local  arrangement  of  molecules  within 
the  monolayer  is  hexagonal  and  the  density  is  higher  than  in  the  bulk  phase.  Good  agree¬ 
ment  was  found  between  the  simulation  results  and  X-ray  diffraction  and  reflection  data. 


Conformational  Analysis  of  Complexes  of  Phosphine  Oxide  Calix-[4]-arenes 
with  Group  1  and  Group  2  Cations 

Paddy  Kane,  Kevin  Kincaid,  Darren  Fayne  and  Dermot  Diamond 

School  of  Chemical  Sciences,  Dublin  City  University,  Glasnevin,  Dublin  9. 

Monte  Carlo  conformational  searches  in  combination  with  energy  minimisation  by  molecu¬ 
lar  mechanics  were  used  to  model  the  complexes  formed  between  two  calixarenes,  1  and  2, 
and  the  Na+,  Li+,  K+,  Mg2+  and  Ca2+.  1H  NMR  titrations  were  carried  out  to  investi¬ 
gate  the  stoichiometry  of  the  complexes  and  the  structural  changes  which  take  place  upon 
complexation.  The  mode  of  complexation  for  each  complex  is  examined  and  the  results  are 
compared  with  the  results  of  selectivity  studies.  2  is  known  from  potentiometric  studies  to 
be  selective  for  Ca2-f  in  the  presence  of  Li-f ,  Na+  and  Mg2+.  1H  NMR  studies  indicated 
that  the  2:Ca2-)-  complex  has  a  1:1  stoichiometry.  The  modelling  studies  indicated  that 
the  Ca2+  ion  was  co-ordinated  by  the  eight  oxygen  atoms  of  2  and  that  this  mode  of  com¬ 
plexation  was  unique  for  the  2  complexes.  In  contrast,  an  eight-fold  mode  of  complexation 
was  observed  in  the  models  of  the  complexes  formed  by  1  with  three  of  the  ions  studied  — 
Na+,  K+  and  Ca2+.  These  results  clearly  demonstrate  why  2  shows  a  strong  preference 
for  calcium  in  the  presence  of  the  other  interferents  studied  and  also  why  1  has  a  poor 
ability  to  discriminate  between  calcium  and  some  of  the  other  interferents. 


GW  Calculations  in  a  Gaussian  Orbital  Basis 
C.H.  Patterson 

Department  of  Physics,  Trinity  College,  Dublin. 

Hartree-Fock  and  Density  Functional  Theory  methods  applied  to  electrons  in  solids  are 
well-known  to  result  in  band  structures  and  other  properties  which  differ  from  each  other 
and  from  experiment.  The  usual  problem  is  that  the  Hartree-Fock  approximation  results 
in  band  gaps  that  are  too  large  and  density  functional  theory  results  in  gaps  that  are 
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too  small,  or  non-existent.  Hedin’s  GW  method  has  been  applied  to  semiconductors, 
insulators  and  metals  in  the  last  12  years  or  so  and  generally  gives  a  satisfactory  account 
of  the  electronic  structures  of  these  materials.  GW  calculations  in  a  Gaussian  orbital  basis 
are  discussed  as  well  as  GW  results  for  an  insulator  in  the  extreme  tight-binding  limit. 


Infrared  absorption  as  a  probe  of  strain  and  microstructure  in  boron  nitride 
thin  films 

Stephen  Fahy 

Department  of  Physics,  University  College,  Cork. 

Experimental  and  theoretical  investigations  of  the  infrared  active,  polarization  dependent 
phonon  frequencies  of  cubic  boron  nitride  films  have  been  performed  [1]  in  light  of  recent 
claims  that  large  frequency  shifts  during  initial  nucleation  are  the  result  of  strain  caused  by 
highly  non-equilibrium  growth  conditions.  We  show  that  the  formation  of  small,  separate 
grains  of  cubic  boron  nitride  during  the  initial  growth  leads  to  a  frequency  shift  in  the 
infrared  active  transverse  optic  mode  polarized  normal  to  the  substrate  which  is  opposite 
in  sign  and  twice  the  magnitude  of  the  shift  for  modes  polarized  parallel  to  the  substrate.  In 
contrast,  film  strain  causes  a  frequency  shift  in  the  mode  polarized  normal  to  the  substrate 
which  is  much  smaller  in  magnitude  than  the  frequency  shift  for  modes  polarized  parallel 
to  the  substrate.  Normal  and  off-normal  incidence  absorption  measurements,  performed  at 
different  stages  of  nucleation  and  growth,  show  that  large  frequency  shifts  in  the  transverse 
optic  phonon  modes  during  the  initial  stage  of  growth  are  not  compatible  with  the  expected 
effects  of  strain,  but  are  in  large  part  due  to  nucleation  of  small  isolated  cubic  BN  grains 
which  coalesce  to  form  a  uniform  layer.  Numerical  results  from  a  simple  model  of  island 
nucleation  and  growth  are  in  good  agreement  with  experimental  results. 

[1]  S.  Fahy,  C.  A.  Taylor,  and  R.  Clarke,  Phys.  Rev.  B  56,  12573-80  (1997). 


Atomistic  simulation  and  ab  initio  computations  for  microelectronics 
Jim  Greer 

Materials  Modelling  Group,  National  Microelectronics  Research  Centre,  UCC,  Cork. 

In  1997,  the  NMRC,  Ireland  established  a  new  unit  within  the  Advanced  Materials  & 
Technology  Group  to  apply  material  computation  methods  to  problems  of  relevance  to 
microelectronic  issues.  Material  computation  work  represents  a  clear  departure  from  the 
TCAD  used  within  the  microelectronics  industry.  The  emerging  use  of  ab  initio  and  atom¬ 
istic  simulations  for  microelectronic  applications  reflects  the  fact  that  lateral  dimensions 
of  a  few  hundred  nanometers  and  vertical  dimensions  of  hundreds  of  angstroms  are  now 
routinely  being  manufactured.  NMRC  material  computations  have  concentrated  on  ap¬ 
plication  of  molecular  dynamics  simulations  to  studies  of  liquid  crystal  polymers  (see  the 
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abstract  by  Mick  Nolan  in  these  proceedings),  investigation  of  potential  barriers  at  strained 
Si/SiGe  interfaces,  calculation  of  defects  in  compound  semiconductors,  software  and  appli¬ 
cation  work  on  a  Monte  Carlo  Cl  code,  and  newly  we  are  moving  into  the  area  of  quantum 
device  simulation  in  collaboration  with  Physics,  UCC.  In  this  presentation,  we  will  discuss 
the  areas  of  application  for  ab  initio  and  atomistic  simulation  in  the  microelectronics  sim¬ 
ulation  hierarchy.  Discusson  of  recent  results  obtained  at  the  NMRC,  and  in  particular, 
calculations  of  the  carbon  tetramer  using  the  MCCI  method  will  be  presented. 


Why  is  subsurface  hydrogen  more  reactive  than  surface  hydrogen  in  catalytic 
hydrogenation? 

A.  Michaelides 

Department  of  Chemistry,  Queen’s  University,  Belfast. 

In  the  catalytic  hydrogenation  of  hydrocarbons,  subsurface  hydrogen  has  been  experi¬ 
mentally  demonstrated  by  Ceyer  and  co-workersto  be  much  more  reactive  than  surface 
hydrogen.  We  use  density  functional  theory  to  identify  low  energy  pathways  for  the  hydro¬ 
genation  of  methyl  adsorbed  on  Ni(lll)  by  surface  and  subsurface  hydrogen.  The  unique 
reactivity  of  subsurface  hydrogen  is  found  to  be  due  to  its  excess  potential  energy  rather 
than,  as  has  previously  been  suggested,  a  more  favourable  transition  state.  The  rate  limit¬ 
ing  step  is  the  resurfacing  of  H.  These  features  are  expected  to  hold  for  other  hydrogenation 
reactions. 


Onset  of  Magnetic  Order  in  Ferromagnetic  Nanostructures 
M.D.  Kuzmin 

Department  of  Physics,  Trinity  College,  Dublin. 

A  simple  phenomenological  model  is  proposed,  which  describes  two-component  nanostruc- 
tural  ferromagnets.  The  entire  space  is  divided  into  two  domains  filled  with  two  distinct 
ferromagnetic  media  whose  bulk  Curie  points  are  Tel  and  Tc2,  Tcl>Tc2.  The  two  ma¬ 
terials  are  assumed  to  grow  into  each  other  at  the  interface  so  that  the  exchange  coupling 
across  the  interface  has  the  same  intensity  as  in  the  bulk.  Importantly,  the  size  of  the  high- 
Tc  domain  is  restricted  to  a  small  value  (dl)  in  either  1,  or  2,  or  3  dimensions.  The  main 
result  is  that  the  ordering  temperature  of  the  system  is  lower  than  Tel.  The  deviation 
from  Tel  in  the  first  non- vanishing  approximation  is  quadratic  in  A/dl  and  independent 
from  Tc2(A?slAisa  universal  exchange  length). 
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Modelling  oxides:  simplifications  from  ab  initio  density  functional  theory 
Mike  Finnis 

Atomistic  Simulation  Group,  Queen’s  University,  Belfast. 

I  point  out  how  the  usual  Hohenberg-Kohn-Sham  total  energy  functional  can  be  refor¬ 
mulated  in  a  way  that  leads  to  the  simplified  models  such  as  the  ionic  model  and  the 
tight-binding  model.  Tony  Paxton  and  I  derived  a  self-consistent  tight-binding  model  [1] 
on  this  basis  which  we  have  applied  to  studying  the  phases  of  ZrC>2. 

[1]  M.W.  Finnis,  A.T.  Paxton,  M.  Methfessel  and  M.  van  Schilfgaarde,  Phys.  Rev.  Letters, 
81,  5149-5152  (1998) 


QMC  for  excited  collective  states 
Mike  Harris 

Department  of  Physics,  UCC,  and  Max-Planck  Institute,  Stuttgart. 

Excitation  spectra  for  the  dipole  and  quadrupole  plasma  oscillations  in  jellium  spheres  are 
obtained  by  the  analytic  continuation  of  imaginary  time  correlation  functions  computed 
by  the  fixed  node  diffusion  Monte  Carlo  method.  The  results  provide  benchmark  data  to 
test  approximate  theories  for  excited  states  (e.g.,  the  random  phase  approximation,  time 
dependent  density  functional  and  the  GW  method).  The  method  is  also  applied  to  the 
calculation  of  the  plasmon  dispersion  relations  for  both  plasmon  and  spin  wave  excitations 
in  the  homogenous  electron  gas  at  various  electronic  densities. 


We  wish  to  thank  the  following  for  their  contribution  to  the  success  of  this  conference: 

United  States  Air  Force  Office  of  Scientific  Research,  European  Office  of  Aerospace  Re¬ 
search  and  Development. 


7 


